Crowding induced entropy-enthalpy compensation in protein association equilibria 
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A statistical mechanical theory is presented to predict the effects of macromolecular crowding 
on protein association equilibria, accounting for both excluded volume and attractive interactions 
between proteins and crowding molecules. Predicted binding free energies are in excellent agreement 
with simulation data over a wide range of crowder sizes and packing fraction. It is shown that 
attractive interactions between proteins and crowding agents counteract the stabilizing effects of 
excluded volume interactions. A critical attraction strength, for which there is no net effect of 
crowding, is almost independent of the crowder packing fraction. 

PACS numbers: 



Protein-protein interactions are important in many es- 
sential biological functions, such as transcription, trans- 
lation, and signal transduction [1] . A lot of progress has 
been made in understanding protein association in dilute 
solution via experiments and simulations [5H5]. Cells, 
on the other hand, contain various macromolecules, e.g., 
DNA, RNA, proteins, organelles, etc., which constitute 
up to 40% of the cell volume [H] . It is thus crucial to re- 
late in vitro experimental or simulation results to those 
in a crowded cellular environment [71413]. 

Several experimental studies have been performed to 
understand protein-protein interactions in a crowded en- 
vironment [T4Tf2i] . Most attention has been paid to the 
steric excluded volume effects of inert crowding agents on 
the formation of protein complexes [25—27] . Very recent 
studies have also started to probe the effects of attrac- 
tive interactions between proteins and crowders on pro- 
tein association [28l431| . These studies have highlighted 
the importance of accounting for cnthalpic effects aris- 
ing from attractive interactions in addition to commonly 
invoked excluded volume effects. It was found that the 
enthalpic effects can actually increase the binding free 
energy (thereby destabilizing the bound complex) in con- 
trast to predictions based on available theoretical models 
that can only capture entropic effects. 

Most theoretical models of crowding are based on 
scaled particle theory (SPT) of hard-sphere fluids [32J 
or its modified versions and have been applied to inter- 
pret experimental and computational results with vary- 
ing success. The failure of these models in several sit- 
uations highlights an important role played by attrac- 
tive crowder-protein interactions. Therefore, a need for 
comprehensive quantitative theory is warranted, that can 
describe the effects of repulsive as well as attractive 
crowder-protein interactions on the protein-association 
equilibria. 

In this paper, we present a theory that can quan- 
titatively predict the effects of macromolecular crowd- 
ing on the protein association equilibria accounting for 
both repulsive and attractive crowder-protein interac- 
tions. The statistical mechanics and thermodynamics 
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FIG. 1: Schematic diagram of the thermodynamic cycle for 
the formation of the Ubq/UIMl complex. The ubiquitin is 
shown in blue while UIM1 is shown in red. 



of a hard-sphere fluid are adapted to yield an approx- 
imate analytical expression for the protein-binding free 
energy in the presence of spherical crowders. Exten- 
sive replica exchange Monte Carlo (REMC) simulations 
have been performed on two distinct protein complexes 
to test this theory. We find that the theory is in excellent 
agreement with simulations over a wide range of crowder 
packing fractions and crowder-protein interactions. The 
theory identifies the region in parameter space (entropy- 
enthalpy compensation line in a two parameter plane) 
separating entropically stabilized area versus enthalpi- 
cally destabilized one. 

Theoretical development. Figure [T] illustrates a ther- 
modynamic cycle that describes a change in the binding 
free energy AF blnd of two proteins due to the presence 
of crowding molecules. This change, AAF bmd , can be 
expressed as the difference in the binding free energy in 
the absence and presence of crowders and is given by 



AAF bmd (<?(>) = AF (</>) — AF (</> = 0) 
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where AF= rowd ((/>), (a G [A,B,AB]) is the solvation free 
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energy of a protein (or complex) a in a crowded solution 
with crowding packing fraction <j>. [For brevity, we will 
omit the superscript "crowd" below.] 

To obtain an expression for AF a (cj)) in Eq. (JlJ for a 
protein or complex a, let U a (r, O) = J2 i£a Ui(n) be the 
overall interaction between a protein a and a crowder, 
where r is the distance between the center of mass of the 
protein and the crowder and £1 the orientational degree of 
freedom, while u» is the interaction between an atom (or 
residue) i of the protein a and the crowder. For a general 
Lennard-Jones(LJ)-type potential for Uj, it is reasonable 
to assume that for given fi, U a (r, fi) exhibits a minimum, 
-e™(fi), at r = r™(fi). Following the Weeks-Chandler- 
Andersen (WCA) theory, we then decompose U a into the 
repulsive and attractive parts as 



rj (r m _ / C/ Q (r,0) + e™(^) r<r2»(n), 

u a ,Te P {r,u) - | Q otherwise, 

77 rr m - / " e ™ ( ^ r < r " (fi) ' 

U "' attl ' J ~ \f7 a (r,Q) otherwise. 
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The solvation free energy, AF a (<p), of the protein in a 
crowded solution can then be divided into two parts as, 



AF a (<f>) = AF Q , rcp (0) + AF aMt ((j,), 



(3) 



where AF, 



a,rcp(att) 



is the contribution from the repulsive 
(attractive) interaction, respectively. 

The repulsive contribution, Af Qirep , is obtained by 
adopting the SPT. The SPT theory provides the free en- 
ergy for solvating a hard-sphere of radius R a in a bath 
of hard-sphere particles of radius R c as, 

AF a>rep = (3y + 3y 2 + y 3 )4> + (4.5y 2 + Zy 3 )4> 2 

+ 3y 3 } 3 - ln(l - <j>), (4) 

where 4> = 0/(1 — 4>) and y — R a /R c . But can we 
represent an anisometric protein with soft-core protein- 
crowder interactions as a hard sphere with an appropriate 
radius R a to capture protein's solvation behavior accu- 
rately? Here we use the Boltzmann criteria to define R a 
as, 
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where the right-hand side represents the volume encom- 
passed by the condition [/ Q . rep (r, f2) > fk^T. Here, we 
use / = 2 that has been used successfully in previous 
studies [33] . 

Using thermodynamic perturbation theory approach, 
the attractive contribution, A.F a a ttj can be expressed as 
(up to the first order), 



AF a , 



(t / a ,att)rc P = / pU a ,a,tt(r,O,)g (r)r drdQ,, 



where p is the crowder number density related to <f> via 
p — 4>/ (4ttR 3 /3), and go(r) is the radial distribution func- 
tion of the hard-sphere crowders between a protein and 
a crowder. Realizing that go(r) has a maximum g™ ax 
at contact and then decays rapidly to unity, we assume 
9o(r) = 5™ ax for r € [r™, r™+A) and 1 for r € [r™+A, oo) 
with A = (2 1 / 6 - l)R c ~ 0.12R C 34 . We then approxi- 
mate Eq.([6]) as, 

AF Q , att » -pe a S a {Sr + (g™ x - 1)A}, (7) 

where e a = (e™)n is the orientational average of e™,S a — 
J[r™(r2)] 2 df2 the surface area around the protein, and 5r 
the attraction range. Note that here we assume Sr > A. 

To enhance the simplicity and practical value of our 
theory, we use the Carnahan-Starling (CS) equation of 
state for a hard sphere fluid to calculate g™ ax . The CS 
equation of state is known to reproduce the thermody- 
namic behavior of hard-sphere fluids from dilute gas to 
near the freezing transition. The CS expression for (7g lax 
is given by 



max 

9o 



(l-0/2)/(l 



(8) 



and only depends on 0. Note that the first term in Eq. Q 
gives a linear order in <f) while the term containing <7g lax 
yields higher order terms. Combining together Eqs. (1), 
(3), (4) ,(7) and (8), one can easily obtain an estimate 
of crowding induced change in the binding free energy. 
Next, we test this theory against REMC simulations of 
two protein complexes in a wide range of crowder sizes, 
packing fractions and interaction strengths. 

Model and simulation details. To proceed further, 
we first introduce our model and simulation procedure 
briefly. A residue-based coarse-grained model is used to 
simulate protein-protein interactions |35j . This trans- 
ferable protein-protein interaction model was shown to 
yield binding affinities and structures for moderate-to- 
weakly interacting protein complexes in accord with ex- 
periments |35[ 136] . Crowding agents are represented by 
spheres interacting via a repulsive potential, u rcp (r) — 

e r ^ . r _2r'+a- J ' wnere ay is the interaction range set 

equal to 6A. The interaction between a residue i and a 
crowder is given by a modified LJ potential as, 

Ui(r) 



4e,. 



(9) 



(6) 



where <7j = Ki + r c with m being the van der Waals radius 
of the residue i. This form of the potential, instead of 
the usual LJ form, ensures that the interaction range is 
independent of r c . If the interaction between a residue 
i and a crowder is repulsive, we take u rop (r) with 2r c 
replaced by cr.;. 

Replica exchange Monte Carlo (REMC) simula- 
tions are performed on two distinct protein com- 
plexes, ubiquitin-UIMl (Ubq/UIMl) and cytochrome c- 
cytochrome c peroxidase (Cc/CcP) following standard 
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FIG. 2: Plot of the overall interaction between the 
Ubq/UIMl complex and a crowder, Uubq/uiMij as a func- 
tion of r for different orientations Q. The solid curves are 
obtained from Eq. |9| with e c = e m i n (£1) and cr, = ro(fi). 



protocols [35]. Repulsive as well as attractive interac- 
tions are considered for the residue-crowder interaction 
with e r = 1.69fceT for the repulsive interaction, and 
e c = 0.15, 0.3, 0.45, 0.GfcB? 1 for attractive interactions. A 
wide range of crowder size r c and "bare" crowder volume 
fraction, (f> = 4:TrN c r^/3V, are considered where N c is 
the number of crowders and V = L? the volume of the 
simulation box of length L. 

Results and Discussion. The spherical crowders in- 
teract with each other via the distance-dependent soft 
repulsive potential given by it rop (r) with a characteris- 
tic size r c and e r = 1.69/cbT. To apply the SPT the- 
ory (4) for calculating the repulsive contribution of the 
binding free energy, it is necessary to obtain an effec- 
tive hard-sphere radius for such crowders. We define the 
effective hard-sphere radius, R C1 of crowders by the con- 
dition, u Tcp (2R c ) = fk^T with the same / as in Eq. (pH). 
This yields R c = r c + ja r where 7 = ^[(^§-) 1/12 - 1J. 
Note that although for e r = 1.69fcBT one has R c ~ r c , 
in general, R c can be different from r c . The effective 
packing fraction <f) is then given by <j> — 4>o(R c /r c ) 3 . 

Figure [2] presents the overall interaction between the 
complex Ubq/UIMl and a crowder at five different ori- 
entations, illustrating a highly anisotropic and asymmet- 
ric nature of the interaction. It shows that the over- 
all protein-crowder interaction follows the LJ shape of 
the residue-crowder interaction, Eq. ([9]), (see the solid 
curves), with a minimum —e m (VL) at r — r m (Q) for a 
given n. However, the longer-distance tails are underes- 
timated by the same formula as evident in the inset. 

The effective radius, i? Q , for a protein a, determined 
by Eq. ([5| depends weakly on r c and e c as shown in Ta- 
ble[Tj For the repulsive protein-crowder interactions, such 
effective radii for proteins and complexes are sufficient 
enough to calculate the change in the binding free en- 



TABLE I: Effective radius, R a , (in A), for the ubiqui- 
tin (Ubq), UIM1 and the Ubq/UIMl complex for r c = 
12,16,20Afor attractive (e c = 0.15, 0.3, 0.45, 0.6fc B T) and re- 
pulsive (rep; e r = 1.69£;bT) interactions 





Ubq 
12 16 20 


UIM1 
12 16 20 


Ubq/UIMl 
12 16 20 


0.15 


14.13 14.39 14.57 


9.82 10.14 10.37 


15.88 16.17 16.38 


0.30 


14.29 14.54 14.72 


9.99 10.31 10.54 


16.03 16.32 16.53 


0.45 


14.36 14.61 14.79 


10.07 10.38 10.62 


16.11 16.39 16.60 


0.60 


14.41 14.65 14.83 


10.12 10.43 10.66 


16.15 16.44 16.64 


rep 


15.18 15.42 15.59 


10.83 11.13 11.35 


16.92 17.20 17.40 



a 3 



b 2 




FIG. 3: Binding free energy, AAF b (4>), for the Ubq/UIMl 
complex as a function of the crowder packing fraction <f>. The 
symbols are simulation data while solid curves are predictions 
from the theory (see the text). 



ergy, AAF bmd , via Eq. Q. Figure [3] shows an excellent 
agreement between simulation results (black squares) and 
the theory (black solid curves) for the Ubq/UIMl com- 
plex for different crowder sizes. As previously reported 
by us and others, the binding free energy decreases with 
increasing packing fraction (f> and decreasing crowder size 
due to the excluded- volume effect. 

Recent studies |2SH3U] have shown that attractive 
protein-crowder interactions can destabilize protein as- 
sociation. Figure [3] shows that indeed as the attraction 
strength, e c , between a residue and a crowder increases 
the binding free energy also increases with the pack- 
ing fraction <j). For example, for a moderate strength 
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e c = Q.Qk^T the change in the binding free energy at 
<t> = 0.3 (close to the physiological condition) is up 
to about 4 U-qT when the protein-crowder interaction 
switches from repulsive (black) to attractive (purple). In 
order to apply our theory, Eqs. Q-Q, to describe the 
simulation data for various e c and r c , we calculate the av- 
erage attraction strength, e a , and the surface area, S a , 
for the individual proteins and the complex. Note that 
e a is proportional to e c while S a is independent of e c . 
Table |Tl] shows these values for different r c . The theory 
predictions are in excellent agreement with the simula- 
tion data in which the attraction range Sr = 5 A(close 
to ay) is used for all the crowder sizes and attraction 
strengths. 

To check whether the theory can be transferable to 
other protein complexes, we calculate the binding free en- 
ergies for the Cc/CcP complex (total 402 residues com- 
pared to 100 residues for the Ubq/UIMl complex) as 
shown in Fig. [4] With the same Sr, the theoretical pre- 
dictions agree remarkably well with the simulation data. 

The data in Figures 3 and 4 show the competition be- 
tween entropic effects of the excluded volume and en- 
thalpic effects by attractive crowder-protein interactions. 
As previously suggested 130] , the enthalpic effects can 
be approximated to be proportional to protein's surface 
areas and our theory here provides its concrete founda- 
tion from microscopic nature of the protein-crowder in- 
teractions. At high attraction strengths, the enthalpic 
penalty for breaking the crowder-protein interactions (at 
the expense of protein-protein interactions) dominates, 
thus increasing the binding free energy. At some critical 
attraction e£ rlt , the two contributions are canceled out, 
and the binding energy in a crowded solution becomes 
equal to that in the absence of crowders (see green trian- 
gles and curve in Fig. 3b). 

It was observed [3D] that the critical attraction, e£ rlt , 
for which the effect of the excluded volume is canceled 
out exactly by that of the attractive contribution, (i.e., 
AAF blnd = 0), is almost independent of the crowder 
packing fraction, cj>. This is owing to the fact that 
AAF blnd is almost linear in <f> for e c considered. To ob- 
tain e£ rlt estimate, we combine Eqs. Q and j7]) and solve 
for e c that satisfies AAF bmd = up to the linear order 



TABLE II: Normalized average attraction strength, e a /e c , 
and the surface area, S a , (in A 3 ) for Ubq, UIM1 and the 
Ubq/UIMl complex 
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UIM1 


Ubq/UIMl 


r c 


e a /e c 








e a /e c 


Sot 


8 


4.56 


6543 


4.01 


4100 


4.61 


7522 


12 


4.71 


9278 


4.07 


6418 


4.75 


10487 


16 


4.79 


12402 


4.07 


9143 


4.85 


13830 


20 


4.85 


15921 


4.04 


12273 


4.91 


17562 




FIG. 4: Binding free energy, AAFt, for the Cc/CcP complex 
as a function of </>. Symbols and curves are same as in Fig. 3. 

in <f). One then obtains, 

£ crit _ AY /AW + 0((t) 2 ), (10) 

where 

AF = 3(y A + y B - Vab) + 3(y A + Vb~ Vab) 

Hva + vI-vab), (ii) 

AW = 3{e A S A + e B S B -e AB SAB)Sr/(4TTR 3 c e c )(12) 

This yields ef'/feT ~ 0.19, 0.27, 0.36 and 0.44 for 
r c = 8, 12, 16, 20A, for the Ubq/UIMl, and 0.28 and 0.35 
for r c = 16 and 20Afor the Cc/CcP, respectively, con- 
sistent with simulation data in Figs. [3] and [4] We can 
also plot e£ rlt as it changes with crowder size r c as shown 
in Figure [5] For crowder-protein attraction values above 
this line, one will observe destabilization of protein asso- 
ciation and stabilization below this line. 

In summary, we have presented a quantitative theory 
for protein association equilibria in a crowded solution 
for both repulsive and attractive crowder-protein inter- 
actions. The theory is based on the statistical mechanics 
and thermodynamics of a hard-sphere fluid. Even though 
proteins are highly anisometric, the repulsive contribu- 
tion to the binding free energy is described well by the 
scaled particle theory of hard spheres. The expression for 
the attractive contribution is obtained by using thermo- 
dynamic perturbation theory and the radial distribution 
function of hard-sphere fluids. The theory is in excel- 
lent agreement with simulation results for the Ubq/UIMl 
and Cc/CcP complexes over a wide range of the crowder 
sizes, packing fractions and attraction strengths. Our 
work provides theoretical foundation of understanding 
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FIG. 5: Enthalpy-entropy compensation lines (i.e., 
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the protein-protein interactions in a cellular environment 
in which proteins and crowding macromolecules exhibit 
non-specific interactions in addition to the excluded vol- 
ume effects. 
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